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ABSTRACT 


Rytov's  method  (method  of  smooth  perturbations)  has  been  used 
extensively  to  solve  problems  involving  wave  propagation  in  weakly 
inhomogeneous  media.  Critics  of  the  method  have  questioned  the 
closure  procedure  employed  in  the  Rytov  theory  and  have  expressed 
doubt  as  to  the  range  of  validity  of  the  solution.  Some  of  the  criticisms 
have  claimed  the  method  fails  because  it  does  not  treat  multiple  scat¬ 
tering  and  that  it  is  limited  in  applicability  to  precisely  the  same  region 
as  the  Born  single  scattering  approximation.  The  criticisms  are  shown 
herein  to  be  at  fault.  A  multiple  scattering  interpretation  of  the  Rytov 
solution  is  given,  thus  lending  support  to  those  proponents  of  Rytov  who 
claimed  an  extended  range  of  applicability  relative  to  the  Born  approxi¬ 
mation. 
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MULTIPLE  SCATTERING  AND  THE  METHOD  OF  RYTOV 


The  classical  approach  to  the  solution  of  problems  involving  wave 

propagation  in  weakly  inhomogeneous  random  media  has  been  via  pertur- 
1  -3 

bation  methods.  In  general,  only  a  first-order  solution  (Born  or  single 

scattering  approximation)  has  been  tractable  and  useful.  This  solution  is 

necessarily  limited  to  propagation  over  short  distances  by  the  require- 

4  5 

ment  that  perturbations  be  small.  To  overcome  this  limitation,  Rytov  ’ 

rewrote  the  inhomogeneous  wave  equation  in  the  field  variable  in  terms  of 

the  logarithm  of  the  variable.  He  then  obtained  the  usual  inhomogeneous 

6  7 

wave  equation  in  the  new  variable  by  a  questionable  closure  procedure.  ' 
Results  obtained  by  Rytov’ s  method  have  been  widely  used  by  workers  in 
radio  and  optical  propagation.  Unfortunately,  however,  Rytov’s  closure 
procedure  appears  never  to  have  been  rigorously  justified. 

g 

DeWolf  has  shown  that  the  Rytov  approximation  is  equivalent  to  the 
Born  and  geometrical  optics  solutions  in  regions  where  the  latter  methods 
are  applicable.  He  suggests  a  far  greater  range  of  applicability  including 

9 

regions  where  multiple  scattering  effects  are  important.  Brown  and 
7  10 

Taylor,  ‘  on  the  other  hand,  have  taken  the  position  that  the  Rytov 
approximation  is  valid  only  in  regions  where  the  Born  approximation  holds. 
Brown,  ^  in  addition,  refers  to  the  Rytov  solution  as  a  single  scattering 
approximation.  It  is  the  purpose  of  this  note  to  point  out  that  (1)  in  a  sense, 
the  Rytov  method  does  treat  multiple  scattering,  and  (2)  on  a  proper  inter¬ 
pretation,  the  Rytov  solution  is  not  restricted  as  claimed  by  Brown  and 
Taylor. 
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Aa  a  starting  point  for  the  discussion,  consider  the  scalar  Helmholtz 


equation  in  one  of  the  field  variables 

[V2  +  k2n2(rj]  A(r)  =  0  , 

where  A(rJ  is  the  field  variable,  is  the  wave  number,  and  n(r)  is  the 
randomly  varying  index  of  refraction.  Let  n{r)  =  1  +  N(rK  Then  we  find 

(V2  +  k2)A(r)  =  -2k2N(r)A(r)  -  k2N2(r)A(r)  .  (1) 

We  now  let 

OD 

A(r)  =  ^  Aj(r)  ,  (?.) 

j=0 

where  Aq(£)  is  the  solution  of 

(V2  +  k2)AQ(r)  =  0  .  (3) 

Combining  (1),  (2),  and  (3)  yields 

00  00 

(V2  +  k2)^Aj(r)=  C-2k2N(r)  -  k2N2(r)]^  Aj(r)  .  (4) 

*  j=l  n=0 

Let 

(T2  +  k2)Aj(r)  =  -2kjN(r)AQ(r)  .  (5) 


The  first-order  perturbation  equation  (5)  is  recognized  as  the  well-known 
inhomogeneous  wave  equation  whose  solution  (Born  or  single  scattering 
approximation)  is 


Ajtr)  =  2kp/dv'G(r,r')N(r')A0(r')  ,  (6) 

v  v 

where  G(r,  r_')  is  the  Green's  function  for  the  operator  (V2  +  k^).  To  obtain 
the  higher -order  scattering  terms,  we  subtract  (5)  from  (4)  to  obtain 

00  00  00 

(V?  +  k0>2  Aj<£.)  =  -2kgN(r)  £  A^r)  -  k2N2(r)  £  A.(r) 
j =2  j=l  j=0 

Now  for  n2  2,  let 

(V2  +  k2)A.(r)  =  -2k2N(r)A._1(r)  -  k2N2(r)A._2(r )  (7) 

Since,  for  the  solution  of  interest 

A.  (r)/A  _2<L)>>n(I> 

J  J 

the  second  term  on  the  right  of  (7)  may  be  neglected.  With  this  simplifi¬ 
cation  and  the  definition 


p  (r)  =  A  (r)/A0(r) 

J  J 

the  solution  of  (7)  may  be  written  as 
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(8) 


^4(r)  7  r  Aft(r') 

pj(-’  =  A^iJ  *  2koiv  dv'  0(r.r')N(r')  Pj.jf) 

On  combining  (3),  (5),  and  (8),  we  obtain  the  complete  perturbation  solution, 
placing  multiple  scattering  terms  in  evidence: 

00 

A(r)  =  AQ(r)  pj(r) 
j=0 

=  A0(£.)  C1  +  Pi<L)  +  (£.)+•••  3  •  (9) 

It  proves  interesting  to  write  similarly  the  Rytov  solution  as  a  power 
series  in  orders  of  the  refractivity  N  and  to  make  a  comparison  with  (9). 

The  Rytov  solution  applies  to  the  equation 

vS(r)  +  [v*(r)]  2  +  k2[l  +  N(r)]2  =  0  (10) 

4r  / \ 

obtained  by  setting  A(rJ  =  eT  —  in  (1).  Rytov  set  1f(r)  =  (»q(£.)  +  (fj(r)  and 
defined  (^(r)  as  the  solution  of  (10)  for  N(:r)  =  0.  He  then  obtained 

V2  (f  j  (r )  +  V^fr)-  [2Vi)0(r)  +  Vi|i1(r_)]+  2k2N(r)  +  ^N^r)  =  0 

and  effected  closure  by  neglecting  the  terms  in  [vtj(r)]2  and  N^r).  The 

5 

resulting  solution  for  fj(r)  is  given  by  (8)  with  j  =  1.  Thus,  we  obtain 


A(r)  =  exp[\|(0(r)  +  (fj(r)]  =  AQ(r)  e*p[pj(r)] 

=  AQ(r )  [1  +  P  j  (£.)  +  P  !<£.)/ 2  +  •  •  •  pjtrj/n!  +  .  . .  ] 

=  A0 (r)  [1  +  £  P-'(r)]  .  (11) 

n=  1  J 

where 

2  r  A0<z‘>  PJr‘y 

Pj'(r_)  =  2k0J  dv'  G(r ,  r 1  )N(r 1 ) ~a~(F]  ~^j! -  *  (12) 

A  multiple  scattering  interpretation  of  the  Rytov  solution  is  evident  upon 
comparison  of  (11)  and  (12)  with  (8)  and  (9).  It  is  noted  that  zero  and  first- 
order  scattering  terms  of  the  two  series  are  identical.  The  second-order 
scattering  term  in  the  Rytov  method,  in  the  context  of  the  successive  per¬ 
turbation  analysis.  is  obtained  by  use  of  the  approximation  that  the  first- 
order  scattered  field  everywhere  in  the  scattering  region  is  given  by  the 
constant  p^(r_)/2,  which  is  half  of  the  first-order  scattered  field  at  the 
observation  point  £.  This  is  evidently  a  reasonable  constant  approximation 
to  the  scattered  field  since  fluctuations  due  to  scattering  increase  mono- 
tonically  with  penetration  of  the  medium.  Similarly,  the  jth-order  scatter¬ 
ing  term  is  obtained  by  use  of  the  approximation  that  the  (j  -  l)th  scattered 
field  is  uniform  in  the  scattering  region  and  given  by  pj  ^ (r_)/ j!  • 

The  comparison  of  (9)  and  (11)  reveals  that  the  Rytov  solution  is  valid 
everywhere  in  the  region  of  validity  of  the  Born  approximation  [p(rj  «  13- 

More  than  that,  it  suggests  a  wider  range  of  validity  as  was  claimed  by 

4  5  8 

Chernov  and  Tatarski.  ’ 
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7 

Brown  has  calculated  a  measure  of  the  error  in  fj(r)  occasioned  by 

Rytov's  closure  procedure.  He  concludes  that  the  error  is  small  only  when 

12 

the  Born  approximation  is  valid.  However,  as  Fried  has  pointed  out. 
Brown  treats  the  solution  for  t  ^(r.)  *n  (11)  as  complete  whereas  it  is  more 
properly  only  a  fluctuation.  The  latter  interpretation  is  necessary  to  satisfy 
conservation  of  energy.  ^  Thus,  we  find 

<e*pHj(r)  +  ♦j(r)J)  =  <exp[2x(r)]> 

=  pcp{2[<x2(r)>  -  <X(r)>2  +  <x(r)>]} 

=  exp[2(oZ  +  m)] 

=  1  ,  (13) 

2 

where  ^s(r)  =  Re{^(£)},  o  i6  the  variance,  and  m  is  the  mean  of  x(r).  We 

15 

have  used  the  result  that  y(£.)  is  a  Gaussian  random  variable  in  the  develop¬ 
ment  of  (16).  But  it  is  apparent  from  (6)  that  (x(£))  =  ®  and  that  (x^(£))  * 
Hence,  to  satisfy  (13),  we  interpret  ^(r)  as  a  fluctuation  about  a  mean  value 
-<X  (£.))•  it  is  to  be  noted  that  the  solution  of  (10)  contains  an  arbitrary 
constant.  In  establishing  the  energy  balance,  we  are  merely  applying  the 
appropriate  boundary  condition.  Thus,  the  proper  form  of  the  Rytov 
solution  is 

A(r)  =  AQ(r)  exp  [-(x2(r)>  +  ^(r)]  .  (14) 
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The  multiple  scattering  interpretation  made  in  connection  with  (11)  applies 

2 

to  (14),  but  now  each  scattered  term  must  be  reduced  by  exp[-(x  (r))]. 
Brown's  conclusions,  which  are  appropriate lo  (11),  cannot  be  said  to  apply 
to  (14). 

Taylor  notes  that  (9)  satisfies  energy  conservation  to  O(N^),  i.e., 

when 


A(r )  =  AQ(r)[l  +  pjfc)  +  <P2(r)>]  ,  (15) 

and  that  a  similar  result  is  not  obtained  with  the  Rytov  solution  as  given  by 
(11).  He  then  concludes  that  (ij(r)  must  be  small  if  the  energy  balance  is  to 
be  approximately  maintained.  Interestingly,  however,  when  the  terms  of 
the  multiple  scattering  series  representation  of  (14)  corresponding  to  those 
in  (15)  are  used,  energy  is  also  conserved.  Thus,  we  find 

[1  -  2<X2(r))]<|l  +  *,«£.)  +  *2(r)/z|2> 

»  [1  -  2<X2(r)>]  [1  +  <|*,<r)|2>  +  <*2(r)+  »*2(r))/2] 

=  [1  -  2<xZ(r»][l  +  2<x2(r)>] 

*  1 

4 

if  terms  of  0(N  )  and  higher  are  neglected. 
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Since  the  arguments  of  .Brown  and  Taylor  do  not  apply  to  (14),  it  does 
not  follow  that  the  range  of  validity  of  the  Rytov  solution  is  that  of  the  Born 
approximation.  The  claims  of  a  greater  range  of  validity  made  by  the  early 
proponents  of  Rytov  seem  reasonable  and,  indeed,  are  supported  by  recent 
experimental  evidence.  ^  The  precise  limitations  remain  inadequately 
explored. 
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Rytov's  method  (method  of  smooth  perturbations)  has  been  used  extensively 
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does  not  treat  multiple  scattering  and  that  it  is  limited  in  applicability  to 
precisely  the  same  region  as  the  Born  single  scattering  approximation.  The 
criticisms  are  shown  herein  to  be  at  fault.  A  multiple  scattering  interpretation 
of  the  Rytov  solution  is  given,  thus  lending  support  to  those  proponents  of  Rytov 
who  claimed  an  extended  range  of  applicability  relative  to  the  Born  approxi¬ 
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